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Abstract 

The main argument against the reality of the electromagnetic vacuum 
fluctuations is that they do not activate photon detectors. In order to meet 
this objection we propose a classical model of a photodetector which, in 
the simple case of a light signal with constant intensity, gives a counting 
rate which is a non-linear function of the intensity. For sufficiently large 
signal intensity,the counting rate is proportional to the intensity, in agree- 
ment with the standard quantum result, but there is a dark rate when the 
signal intensity is low. 

PACS numbers: 03.65.Bz, 42.50.Dv, 42.50.Lc. 

1 Introduction 

The aim of this paper is to present a model of a photon counter resting 
upon classical electromagnetic (Maxwell's) theory, but assuming that the 
whole space is filled with a random radiation field (zeropoint field, ZPF). 
The ZPF is assumed to have an average energy ^hu per normal mode, u 
being the frequency and h Planck's constant. Thus the ZPF corresponds 
precisely to the vacuum fluctuations of quantum electrodynamics, taken 
as a real field, and light signals would be radiation above that "sea" of 
ZPF. The reader may ask what is the use of the model which we propose 
because, according to the standard wisdom, classical Maxwell's theory has 
been superseded by quantum electrodynamics and there exists already a 
well established quantum theory of photon detection (see, e.g. the book 
of MandelpP). In order to answer the question the present section will be 
attempt to explain the relevance of our model. 

More than eighty years have elapsed since the discovery of quantum 
mechanics, but a vivid discussion still exists about its interpretation. The 
various interpretations may be roughly classified in two groups which 
could be called realistic and pragmatic. According to the pragmatic view, 
physics just consists of a set of rules for the prediction of the results of 
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experiments. It is assumed that the question whether these rules lead to 
a picture of the world does not belong to the domain of physics but rather 
to the realm of philosophy. This view was supported by Niels Bohr and 
gave rise to the "Copenhagen interpretation", which has been accepted 
by the mainstream of the scientific community until recently. In contrast, 
the realistic view claims that the essential aim of natural sciences, physics 
in particular, is to provide knowledge about the material world, the pos- 
sibility of predicting the results of experiments being a by-product. Of 
course it is an important by-product, because on it rests the criterion for 
the validity of the theory, that is the agreement of the predictions with 
the empirical results. The realistic view was maintained by several of the 
founding fathers, in particular Albert Einstein, and we strongly support 
it. Actually the realistic interpretation follows the tradition of all sciences, 
not only physics, until the advent of quantum mechanics. If the Copen- 
hagen interpretation has had such a big success it is because people have 
failed to find a general and coherent realistic interpretation of quantum 
mechanics. 

The difficulty for a realistic interpretation of quantum mechanics is 
that it seems incompatible with the quantum formalism, at least if real- 
ism is combined with the demand of locality, an incompatibility proved 
by John Bell more than forty years agol2]. During the eighty years of 
life of quantum mechanics very many of its predictions have been tested 
and the agreement between theory and experiments has been spectacu- 
lar, especially in quantum electrodynamics. Thus it is not strange that 
most people are so fond of quantum mechanics that Bell's theorem has 
been taken as a proof of the impossibility of a local realistic theory of the 
world, even without the need of any additional empirical support. This 
prejudice has produced a bias in the analysis of performed experiments, 
leading to the current wisdom that local realism has been empirically re- 
futed, modulo a few irrelevant loopholes[4]. However, the truth is that 
all experiments performed till now have had results compatible with lo- 
cal realism and the alleged empirical disproof rests upon some unjustified 
extrapolations [3]- Of course quantum mechanics has not been violated, 
so that local realism and quantum mechanics are compatible for all ex- 
periments performed till now. Therefore it is an open question whether 
they are complatible for all experiments which can be actually performed. 
If this is true, then some of the assumptions needed to prove Bell's theo- 
rem are flawed. This means that either the quantum formalism or, most 
probably, the quantum measurement theory should be modified. 

The goal of devising a modification or alternative to quantum mechan- 
ics compatible with local realism (and the experimients) is not easy. A 
possible procedure is to begin with some theory which: a) provides from 
the start a clear picture of the natural world, and b) allows an interpreta- 
tion of the experiments within some restricted domain. The said theory 
should be generalized later in order to cover a wider domain. A theory of 
this type is stochastic electrodynamics (SED)^ whose basic idea is that 
quantum behaviour is related to the existence of a real ZPF. 

The existence of vacuum fluctuations is a straightforward consequence 
of field quantization [6]. In addition, quantum vacuum fluctuations have 
consequences which have been tested empirically. For instance, the vac- 
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uum fluctuations of the electromagnetic fleld give rise to the main part of 
the Lamb shift [3 and to the Casimir effect [S]. Thus Planck's radiation 
law should be written in the form 



huj 1 ^ 



exp{huj/kBT) - 1 2 



^coth— , (1) 



where the second term represents the ZPF. That the thermal spectrum 
contains an lu^ term has been proved by experiments measuring current 
fluctuations in circuits with inductance at low temperature Of course, 
the term is ultraviolet divergent so that some cutoff should be assumed, 
possibly at about the Compton wavelength, where the fluctuations of the 
Dirac electron-positron sea become important. The ZPF of SED is iden- 
tical with these quantum vacuum fluctuations, taken as a real stochastic 
fleld. 

The standard wisdom is that the vacuum fluctuations cannot be inter- 
preted as a real random electromagnetic fleld because they do not acti- 
vate photodetectors in the absence of signals. There is also a gravitational 
problem because, if the quantum vacuum fluctuations are at the origin of 
the cosmological constant, as is usually assumed, that constant should be 
many orders of magnitude larger than the observed value, but we shall 
not be concerned with gravitational effects in this paper. A common 
explanation of the fact that the zeropoint fleld (ZPF) does not activate 
photodetectors is to say that the ZPF is not real, but virtual. In our 
opinion replacing a word, real, by another one, virtual, with a less clear 
meaning is not a good solution. In the present article we shall prove that 
the behaviour of photodetectors can be explained without renouncing the 
reality of the ZPF. The proof goes via constructing an explicit model a 
of detector, producing a counting rate proportional to the intensity of the 
signal, that is also able to subtract efficiently the ZPF. The mere exis- 
tence of the model proves that the said objection to the ZPF reality is 
untenable. Such a proof is the first purpose of the present paper. 

In addition we think that a classical model of detection may add to 
understand the origin of an open problem lasting for more than 30 years, 
namely the existence of a detection loophole. In fact, as is well known 
all experiments aimed at a discrimination between quantum mechanics 
and local realism, via the test of Bell's inequalities using optical photons, 
suffer from a loophole due to the lack of efficient detectors. The standard 
wisdom is that the difficulty in manufacturing efficient photon counters 
is a minor technical problem. However the persistence of the difficulty 
for so many years may indicate that the problem is a fundamental one. 
The comparison of the classical model that we propose with the quantum 
detection theory may throw light on the problem, and this is another 
purpose of the present paper. 



2 Stochastic properties of the zeropoint 
field 

Now we shall derive some relevant properties of the ZPF in free space, that 
is far from any material body. The ZPF is characterized by the electric 
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field E(r,t), and the magnetic field B(r,t). The stochastic properties of 
the field may be summarized saying that it is Gaussian, so that the mean 
and the correlation functions are sufficient to characterize all stochastic 
properties. The mean is zero and the correlation functions of the field 
components arejlOj pTj 

{Mr,t)EAr',t')) = I d'k (^S., - cos [k-{r-r')-c.{t-t')], 

{B,{v,t)B,{v',t')) = {E,{v,t)E,ir',t')), 

{E,Xv,t)B,{v',t')) = |d»te,^,||^cos[k.(r-r')-^(t-i')]. (2) 

where k = |k| = cj/c. The integrals in these expressions do not converge 
but a natural cut-off in u appears as explained in the following. Our aim 
is to make a photodetection model and we may assume that any detection 
event takes place essentially at a single atom. Thus the quantity of interest 
for us will be the time autocorrelation of the averages of the fields over 
the volume of the atom. We define the averaged quantities 

E, (t) = I d\E, {r,t)p (r) , B, (t) = J d\B, (r, t)p (r) , (3) 

where no confusion should arise by the use of the same label for the field 
at a point and the averaged quantity. The function p (r) represents a nor- 
malized effective electron density in the atom, which I assume spherically 
symmetrical. From eqs.ll^} it is straightforward to get the time correla- 
tions of the averaged quantities and we get 

{E,{t)Ek{t')) ^ {B,{t)Bu{t')) ^ 5,kF{t-t'), {E,{t)B,it'))=0, (4) 

where 



32-Kh ^ 



F(t) = / uj cos (lot) duj 

3c 



r 

Jo 



sin {lot /c) p (r) rdr 



(5) 



For any reasonable density, p{r), the function F{t) — F (—t) decreases 
slowly for ]r] < a/c and rapidly for higher t, becoming negligible for 
|r| >> a/c, c being the velocity of light and a the atomic radius. 

For our detection model the most relevant quantity is the radiation 
intensity represented by the Poynting vector, 

S{t) = ^E(t) X B(t), 

where E(t) and B(t) are the averages over the volume of the atom de- 
fined in |(3]). The mean of the Poynting vector is zero and the correlation 
functions of its components may be derived from @ to be 

{S,{t)S,{t')) = ^5,,F {t' - tf ^ a'S,,S {t' - t) , (6) 

where Sjk is Kronecker's delta, 5 (t' — t) is Dirac's delta and o" is a con- 
stant of the order of 10~^c^''^a~^'^^. The second equality involves approx- 
imating every component of the Poynting vector by a white noise, which 
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is plausible provided we study detection of light signals whose coherence 
time, T, is much larger than the time interval where F{t)^ is not negli- 
gible, that is r >> a/c. These properties characterize the components of 
the Poynting vector as three stationary stochastic processes, statistically 
independent and having the autocorrelation of a white noise. 

Actually we are interested in the ZPF within a photon detector and the 
question is whether we may use there the statistical properties ([BJ of free 
space. Our answer is in the affirmative because we should assume that the 
ZPF permeates everything. Indeed in SED it is supposed that the vacuum 
fluctuations of the electromagnetic field - and possibly other fields also - 
is the cause of the random position of the electron inside the atom. In 
the standard interpretation of quantum mechanics it is also assumed that 
the vacuum fields fluctuate even inside material bodies, with an average 
energy ^hu per normal mode, although the modes are modifled by the 
presence of matter (this gives rise, for instance, to the Casimir effect and 
the Lamb shift as mentioned above) . However the radiation modes which 
are most modified by the presence of matter are those of low frequency, 
much lower than the ones most strongly involved in F{t), eq.Q . Thus 
we may use the free space statistical properties ((6} in our model. 

We shall study the situation where we have a determinate light sig- 
nal superimposed on the ZPF. The statistical independence of signal and 
random fields leads to (compare with ([61 ) 



{S,(t)) = Sjals, {Sjit)Skit')) = S.sSksIs + ry^SjkS {t' - t) , (7) 



where S3 is the component of the Poynting vector in the direction of 
the beam and Is the signal intensity. Eq.Q will be the basis of our 
subsequent study but we should bear in mind that the approximations 
involved may be too crude for some applications. Possible improvements 
will be considered elsewhere. 

3 Detection model 

Several models of photon counter have been proposed in the context of 
SED, resting upon the idea that there exists a "detection time", T, in- 
dependent of the light intensity and such that the probability of a count 
depends on the radiation, including the ZPF, which enters the detector 
during the time T[T2]. Those models, however, are not compatible with 
empirical evidence|13). Instead of fixing the detection time T, here we 
shall assume that every atom accumulates the energy and momentum of 
the radiation arriving at it and a count is produced whenever that energy 
surpasses some threshold. If S(t) is the Poynting vector of the radiation 
arriving at the detecting atom at time t, the accumulated energy at time 



where A is the effective cross section of the atom. In the following we 
shall put yl = 1 in calculating the detection rate and multiply times A at 
the end. 



T will be 
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The essential assumption of our model is that a detection event is 
produced at a time T, after the previous count, when T is such that 

E{T) = |E(r)| = E^, E{t) = f S{t')dt' , (9) 

Jo 

and E is a parameter characteristic of the detector. 

The calculations resting upon the model Q are involved due to the 
fluctuations of the ZPF and the signal. Indeed constructing a detailed de- 
tection model on the basis of that equation would require using the theory 
of "first passage time" for the (vector) stochastic process E(t), which has 
a finite, nonzero, correlation time. The problem is dramatically simplified 
if we assume that every cartesian component of E(t) is a white noise (hav- 
ing a null correlation time) superimposed on a deterministic signal with 
constant intensity Is, as in eq.([7]). More specifically, we assume that every 
component of the stochastic process E(f) is a Wiener (Brownian motion) 
process. The calculation of the first-passage time is now straightforward. 
We must begin solving the diffusion equation 

dp I 2 / d^P d^P d^p\ T dp N /.,„N 

with the initial condition 

p (E, 0) = 5^ (E) = 5 (£,) & (Ey) S [E,) , (If) 

and an absorbing barrier at jEj = Em, that is the boundary condition 

p(E,t) = Ofor 1E| ^Em. (12) 

The solution of ea. (|10|l with these conditions is cumbersome due to the 
fact that the symmetry properties of the boundary condition, eq. (|12p , are 
different from those of the diffusion eq. ((10} itself. In order to simplify the 
calculations, whithout changing qualitatively the model, we substitute for 
eg. flf 2|l the following 

p (Ex, Ey, Ez,t) = for either E^ — ±Em or Ey = ±Em or Ez — ±Em- 

(13) 

This means replacing a sphere in the 3D space of the vectors E by a cube. 
Thus our problem may be reduced to solving 3 one-dimensional equations 
by introducing the functions fx, fy, fz such that 

p{Ex,Ey,Ez,t) = fx{Ex,t) fy[Ey,t) fz[Ez,t) . 

Indeed the solution of ea. (llO|l with conditions eas. (|ll|l and p3|) may be 
obtained from the solution of the three ordinary differential equations 

df^^l 2^ df^^l dh^l 2d^_j df^ 

dt 2"" dt 2*^ dEy dt 2" dEi dEz' ^ ' 

with initial and boundary conditions, respectively, 

fx {Ex,0) =S{Ex), fx {Em,t) = fx {~Em,t) = 0, 
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and similar for fy and /z. The appropriate solutions of eas. (|14|l may be 
obtained by the method of the images and we get (we shall put cr^ = 2 in 
the calculation and restore a in the final result) 



{E, + 2nE^f _ Ilt_ 
U 4 



(15) 

and similar expressions, except replacing Is by 0, for fx {Ex,t) and fy(Ey,t). 

Now our model predicts that, if we have a detection event at time 
i = 0, the probability that the next detection event takes place before 
time t is 

/Em rEm j-Em 

fx{Ex,t)dEx fy {Ey,t) dEy / {E,,t) dE, . 

- Em J-Em J-Em 

(16) 

Our aim is to calculate the detection rate, which is the inverse of the mean 
first passage time, that is 

<t>=rt^rf. (17) 

The proof that this average gives the inverse of the detection rate is as 
follows. We consider that the detector is active during a very long time 
interval. Within it we will have a large number of detection events on 
each atom. Let us assume, for the sake of clarity, that the time intervals 
between two detection events form a discrete sequence ti, t2, ...tj,... If we 



llN.t, Y.Pjh <t>' ^ ' 



events has duration tj. If we pass to the continuous, we shall replace the 
summation by an integral, giving a rate R equal to the inverse of < t >, 
which completes the proof. The counting rate in a macroscopic detector 
consisting of many atoms may be obtained from the predicted rate in one 
atom but this will not be made in the present article. 

According to eas. (|16p and (|17p the mean "first passage time" is given 

by 

< i>=/ dtLK^,L= f,{E,,t)dE, (19) 

Jo J-Em 

/Em r-Em 
fx{Ex,t)dEx= / fy(Ey,t)dEy. 
'Em J-Em 



Hence the calculation is straighforward putting eas. (|15|l in eq. (|19|) and 
performing the integrals. But before proceeding at this calculation we 
will work a one-dimensional model consisting of the replacement of K{t) 
by unity, which corresponds to substituting a layer for the atom. As we 
shall see that model gives predictions not too different from the three- 
dimensional one and it has the advantage of leading to a simple analytical 
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expression for the detection rate. From eas. (|15|l and (|19|) we obtain, mak- 
ing the time integral first, 



< t > = 



dtL 



Em 1 

dE —= 

-Em V4jT 



exp 



:EIs 



(20) 



J n 



dt—p exp 



[E + 2nE„ 
it 



■ exp 



-Is \E + 2nE„ 



where the integral has been taken from Gradshteyn[T5] . In order to 
perform the sum in n it is convenient to separate terms with n = 0, n > 
and n < 0, giving 



Sum 



exp 



— exp 



1 



Is\E\ 



+ E(-ir exp 



-Is{2nEm + E) 



+ exp 



-Is{2nEm-E) 



2(1 + exp [Is Em]) ^cosh 



zLE 



The integral in E is now easy and we obtain 



< t >-- 



E„ 



tanh ( l-IsEn 



and, restoring the parameter a, 



R = 



<t> E„ 



■ coth 



Is Err 



(22) 



(23) 



which has the limits 



R~-^ (1+2 exp -2 

Em 



. Is Err 



„2 „2 2 

for Is » — , R-^-F^ for Is « — . 

(24) 



In order to solve the 3D model we begin by calculating 



K{t) = 



1 



1 



dE, E ("l)"exp 



E 



n— — oo 
(2n + l)B„ 

(2ri-l)E,„ 



[E, + 2nErr 



(-l)"exp 



s 

"it 



4t 



ds, 



(25) 



where the last expression is obtained by means of a trivial change of 
variable. We see that K{t) is the integral of a Gaussian times a function, 
F{s), taking alternately the values +1 and —1 in intervals of length 2iJ,„. 
We will perform the integral using the Fourier expansion of this function, 
that is 

oc . 



{2k + 1) TT 



f{2k + l)TVS 

( — 1) cos ' 



whence we get 



K{t) = 



^V4^t^^2k + l 



{2k + l)ns 



2E„ 



4 (-I)'' 



{2k + ifn^t 

4^2 



fc=0 

(26) 



]} 



(21) 
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After that we have, from eas. (|26p and (|19|l . 

dt 



< t >- 



f 

Jo 



K{tr 







/ dE exp 

J -Em 





E (-i)"exp 



n— — oo 



U 4 



Making the time integral first we get (compare with ea. (|20|l l 

°° °= rEm 1 



■EE 



dE —= 

Em V47r 



exp 



(27) 



where 
Hence 



dt—= exp 



4f 



\2 r2 . 

4 



n — — oo 

■ exp 



(2A: + 1) (21 + 1) 



[(2fc + l)' + (2Z + l)']7r2 



11 + 



£■2 



,,_n ,_n ^ -"fe' ^-Em J 



n— — OO 



1 



-Jfci IS + 2n£;„ 



(29) 



Iki\E + 2nE„ 



The sum in n is similar to that in (|2H) with Iki substituted for la , the 
integrals in E are straightforward and we obtain 



(-1)" 



k+l 



^^SS[^«-^'H2fc + l)(2/ + l) 



^ cosh [\laEm) 

cosh{\lkiEm) 



128 S4 



'-F 



IsE„ 



n-) - EE 



(-1) 



k+l 



G{x) = 1 



,4 [(2fc + ly + (2Z + 1)1 {2k + 1) (2/ + 1) 
cosh X 



G{x), 



cosh Jx^ + iyr^ [(2fc + 1)' + {21 + 1)' 



(30) 



where we have restored the parameter a (see eq.([T5}). The function F{x) 
is defined by an infinite series which should be calculated numerically. 
Here we shall make only the calculation for large and small values of x. 

For 2: >> 1 we may write (using the approximation cosh a ~ i exp a 
, valid for a >> 1) 



G{x) 



1 - exp <J s - yjx'^ + -7r2 [{2k + ly + {21 + 1)' 



1 — exp 



TT^ [(2fc + l)^ + (2/ + l)'] 

82; 



+ 



(31) 



Hence it is easy to get F{x) for large x, that is 



^' 8.^„^(2fc + l)(2; + l) 8x[2^^^{2k + l)) 



128a;' 



Jnkl 1 



(28) 
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whence we obtain 



which is similar to the resuh of the one-dimensional model but the de- 
viation from the standard quantum result decreases as instead of 
exponentially (see ea. (l24|l .) 

For x = we get, from eqs. pO)) and ((29)) , 



OO oo -L y 

no) = EE- 



1-1/ cosh |W(2fc + l)^ + (2/ + l) 



(2fc + 1) (2? + 1) [(2fc + 1)2 + (2; + 1)^] 
The series converges quickly and a numerical calculation gives 

<t >~ 0.49:^ 2.0^, 

that is a rate at zero intensity about twice that in the one-dimensional 
model, ea. p4)l . 

It is remarkable that we obtain a detection rate proportional to the sig- 
nal intensity, that is a perfect subtraction of the ZPF, when the intensity 
of the signal is big, a result in agreement with the quantum mechanical 
prediction. However for low, or nil, intensity signals there is some count- 
ing rate, which we should interpret as a "fundamental dark rate" of the 
detector. 

The result may be generalized to the case where the signal intensity 
is not a constant, but a known function of time. It would be enough to 
substitute Ia{t')dt' for 1st in the above equations, although the solution 
of the differential equation (|10p would be more involved. More difficult 
would be to treat the case where the signal itself fluctuates (with a corre- 
lation time of the order of the inverse of the frequency bandwidth). We 
shall not study these problems here. 

We might also analyze coincidence counts in two detectors when the 
incoming beams, with intensities Iiit) and l2(t) above the ZPF, are cor- 
related. The calculation would be straightforward, although lengthy. We 
may assume that ea. (|24p is still valid for each detector and the coincidence 
rate, with a time delay r, will be 

Ru {h{t)h{t + t)) , (32) 

again in agreement with the quantum prediction. However the current 
situation may not be that. In practice the signals may be stochastic and 
not independent of the ZPF. The calculation in these conditions would be 
rather involved. 



4 Discussion 

Our analysis shows that quantum vacuum fluctuations of the electromag- 
netic field (or ZPF) may be efficiently subtracted by a model which as- 
sumes that the radiation is a classical (Maxwell) field including a fluctuat- 
ing ZPF, provided that the fluctuations of the signal have a large enough 
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correlation time in comparison with the effective correlation time of the 
ZPF. This is usually the case in astronomical observations. In contrast, in 
standard quantum optical experiments the fluctuations of the signal may 
have a rather short correlation time. If the correlation time of the signal 
does not fulfil the assumptions of the previous section, the presence of the 
ZPF will probably give rise to deparatures from the standard quantum 
predictions eas. (|23|l and (|32}, that is they will produce some nonidealities 
in the behaviour of optical photon counters. This is specially important 
when it is necessary to measure coincidence counting rates with short time 
windows, as is frequent in quantum optical experiments (e.g. optical tests 
of Bell's inequality). If this is the case, our approach may provide an 
explanation for the difficulties of performing loophole-free tests of Bell's 
inequality using optical photons. 

I emphasize that, although our model is semiclassical, probably its 
main properties might be reproduced by a rigorous quantum treatment. 
Furthermore the difficulties for reaching an intuitive picture of how detec- 
tors subtract the ZPF probably do not derive from quantum theory itself, 
but from the use of approximations like first-order perturbation theory or 
taking the limit of time t — > oo in calculating the probability of photon 
absorption per unit time. Indeed I conjecture that excesive idealizations 
might be at the origin of the difficulties for undersanding intuitively the 
paradoxical aspects of quantum physics. Although simplifications are ex- 
tremely useful for calculations, they tend to obscure the physics. 



5 Conclusion 

We want to compare our rate prediction 



= = 4^ coth (^-^] , (33) 

<t> Em V Cr / 

with the quantum prediction 

Rg = vKIs, (34) 

where rj is the quantum efficiency and K is a, dimensional constant taking 
into account the energy of one "photon" , the effective cross section of an 
atom and the number of atoms in the detector. The two predictions are 
identical for high intensity provided we identify 

E^. - ^- (35) 

But for low intensity and/or high efficiency there is a departure. For 
instance the departure becomes more than 10% (that is the "fundamental 
dark rate" is more than 10% of the total rate) when IsEm/o-^ < 1.5. Using 
(I35[) this gives 

Is < I.^Kg^t), 

that is when the intensity is too low and/or the quantum efficiency too 
high. As the tests of Bell's inequalities requires both high efficiency and 
low intensity our model may explain the difficulty, or maybe impossibility 
of such tests . A more quantitative result would require an estimate of K, 
which is not easy. 
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6 Appendix 



6.1 Evaluation of F(r) with p (r) = (Stto'^) ^ exp {—r/a) 

The density p (r) is normalized in the sense p (r) A-Kr^dr = 1. 
We start calculating 



^0 



I (uir/c) p (r) rdr — '^-^ ^ — = lix), 



Gw=s/ ^'r:":r ^ (36) 



where we have introduced the dimensionless variable x = ua/c. Thus we 
get ^ 

F(t-t')^^^^ n xcos{Tx)[I{x)fdx, 

Introducing the dimensionless parameter r = \t — t'\c/a and the new 
function 

G(r) = G{\t-t'\c/a) = F{t-t') 

we get 

2hc f °° x^ cos (tx) dx 
W io (x^ + 1)' 
It is possible to get an analytical expression for G (r) (an integration by 
parts may transform it in the integral n° 3.773/1 of Gradshtevn|15)) but 
we shall just obtain the behaviour for large and small r. 
For r << 1, that is \t — t'\ « a/c, we may approximate 

/ \ 1 1 2 2 
cos (rx) ~ 1 — -r x , 

and we obtain, to order r^, (with the change x^ — z) 

he f°° z{l- \t'^z) dz he 2^ 



G{r) 



For r >> 1 the main contribution to the integral will come from 
values of x close to the minimum of the quantity y = x~'^ (x^ + l) , that 
is a; = ^3/5. Making the change of variable x = u + ^3/5 we may 
approximate, neglecting terms of order u^ and higher, 

8* /s"/, , 25 2 



y iu) ~ 7t\ — I 1 H u 

^ ^ ' 3.53 V 3 V 8 

Thus we get, extending the integral in u to -00 because only values of x 
close to ^3/5 will contribute. 



G(i 



5 3 he 



cos 1^ (^u + ^3/5^ r j du 



256 V 5 7ra4 J _^ (u^ + 8/25) 

3" he ( f3^\ cos (tu) du 



256 V5 7ra4 \^V5 ^ (u^ + 8/25) 

25 [Yhe I 2^2 \ ( 

-T COS \ —T 



512 V 10 a* V 5 / \ V 5 
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6.2 Evaluation of a 

From eq.® it is natural to assume 



„2 foa 



o 2 / p^^ydt- 

The integral may be obtained from eq.(|36} performing the r-integration 
first. Taking into account that 

cos(a;r) cos {x't^ dt — -k \5 {x — ;c') + <5 (;c + x')] , 



2 



we get 
giving 



Jo 



3 (47ra 



7/2 - --^^ „7/2 
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